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The low field phase of the organic superconductor (TMTSF)2C1C>4 is studied by muon spin ro- 
tation. The zero temperature limit of the magnetic penetration depth within the TMTSF layers is 
obtained to be A a i>(0) = 0.86(2) /im. Temperature dependence of the muon spin relaxation shows 
no indication of gap nodes on the Fermi surface nor of any spontaneous fields due to time-reversal- 
symmetry breaking. The weight of evidence suggests that the symmetry of this low field phase is 
odd-frequency p-wave singlet; a novel example of odd-frequency pairing in a bulk superconductor. 

PACS numbers: 74.70. Kn, 74.20.Rp, 76.75.+i 



The Bechgaard salts 1 (TMTSF) 2 X (anion X typi- 
cally being PFg or CIO4) attract continuing interest as 
a system whose rich physics is derived from a quasi- 
one-dimensional character and strong electron-electron 
interactions^. Although (TMTSF) 2 C10 4 (TMC) was 
the first ambient pressure organic superconductor (SC) 
to be discovered^, the exact nature of its unconventional 
SC phases has yet to be resolved, with recent studies in- 
dicating that there is more than one distinct SC phase 
in this system. One phase is observed up to poH C 2 = 
0.16 T for magnetic fields H oriented perpendicular to 
the TMTSF layers^ (low-field L-SC phase), but aligning 
H within the TMTSF layers increases H C 2 significantly 
above the Pauli limit 7 and allows a further high-field SC 
phase to be revealed. NMR finds clear evidence^ for a 
field-induced transition between the L-SC phase and the 
second (H-SC) phase for fields above ~1.5 T; the L-SC 
phase shows the Knight shift expected for a singlet state, 
whereas the H-SC does not, indicating a triplet state. 

A fundamental parameter for a superconductor is the 
magnetic penetration depth A, which determines the su- 
perfluid stiffness p — c 2 /A 2 characterising the electro- 
magnetic response to an applied magnetic field. The 
variation of p(T) or A(T) with temperature T provides an 
important test of the gap symmetry. Previous estimates 
of A a b(0) for TMC using muon-spin rotation (/xSR) have 
either given a value around 1.2 /J,m2r— or a distinctly 
smaller value in the region of 0.5 p,m-£. None of these 
studies took into account the strong H variation of the 
line width due to the small value of H C 2 when the su- 
percurrents are in the ab plane. In order to resolve this 
A a b(0) discrepancy and gain further information about 
the SC gap symmetry, we have made a detailed fxSR 
study of this system, investigating the dependence of the 
SC properties over a full range of H and T. 

From these results and other previously reported 
properties, we find strong evidence that the symme- 
try of the L-SC phase in TMC is an unusual odd- 
frequency pairing odd-parity singlet state. Odd fre- 
quency pairing^ 3 - was first proposed^ 4 - for 3 He and subse- 



quently for superconductors—. It is generally believed to 
be most relevant to symmetry-broken local environments, 
such as those near interfaces or in vortex coresi 3 - and has 
not previously been verified in a bulk superconductor. 

Crystals of TMC form as needles aligned with the 
molecular stacking direction a. This is the most con- 
ducting direction of the material, which has a triclinic 
structure in which layers of TMTSF are arranged in the 
ab plane, separated along the c axis by layers of the CIO4 
anions. The c* direction perpendicular to the molecular 
planes is the axis with the smallest conductivity. Our 
sample consisted of a mosaic of these crystals (total mass 
124 mg) aligned with their a axes parallel. In the mea- 
surements H is applied perpendicular to a, so we measure 
an axial average of the properties. 

Transverse field muon spin rotation (TF /iSR) provides 
a means of accurately measuring the internal magnetic 
field distribution caused by the vortex lattice (VL) in a 
type II SCi£ and has previously proven useful in estab- 
lishing the vortex properties of molecular SCs^ 7 -^ 8 . Our 
field-cooled TF /xSR measurements used two /i + SR in- 
struments. At the Swiss muon source (S//S) within the 
Paul Scherrer Institute we used the LTF dilution refriger- 
ator instrument to measure down to 20 mK. At the ISIS 
Facility, STFC Rutherford Appleton Laboratory, UK we 
used a 3 He sorption cryostat with the MuSR instrument 
to measure down to 300 mK. For the LTF measurements, 
magnetic field -B a pp = f^oH was applied perpendicular to 
both a and the plane of the sample. For the MuSR mea- 
surements, -Bapp was perpendicular to a and parallel to 
the plane of the sample mosaic. To ensure that the sam- 
ple was in the non-magnetic fully relaxed anion-ordered 
(AO) state the cooling rate was maintained at 1 K h _1 
to 2 K h _1 between 35 K and 15 K. All data analysis was 
carried out using the WiMDA program 1 ^. 

In TF /iSR measurements performed above T Cl the 
distribution of field B at the muon, p(B), is broad- 
ened mainly by randomly oriented nuclear moments near 
the muon stopping sites; in this case dominated by the 
methyl groups at the ends of the TMTSF molecules. This 
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FIG. 1. (a) SC line width -Brms and (b) average field shift B s h 
versus T for several different values of applied field measured 
using LTF. Characteristic transition temperatures T' are ob- 
tained from two- fluid model fits (lines), (c) The low field 
phase diagram, comparing T' from /iSR (points) with pre- 
viously reported^ T C (H) behaviour for the H\\c* orientation 
that dominates here (dashed line). T' from B Ilns is consistent 
with a 3D vortex melting curve (solid line)2£ 

leads to an essentially Gaussian relaxation of the muon 



polarisation P x (t) oc 

T 2 



cos(7 M (-B)t + 4>), where 

*k = ll(( B - (B)?)- For ^a PP = 2.5 mT the TF 
broadening cr n /j^ was monitored on cooling through the 
AO region and the width was found to reduce by about 
10% on ordering. The AO doubles the b axis leading to 
two inequivalent TMTSF stacks and significant distor- 
tion of the methyl groups on one of the stacks leading to 
a non planar conformation 2 ^. The observed magnitude 
of change in er n is consistent with this distortion. 

In the SC state the VL contributes a further broaden- 
ing, giving a 2 = <r 2 + <jy L with corresponding VL field 
width B lms = ovl/T/j- Due to the relatively weak VL 
contribution, we did not attempt to spectrally resolve 
the asymmetric VL field profile as we did in previous 
studies of organic SC such as ^(ET^CutNCS)^^ (k- 
ET). The T dependence of B lms and the shift of the 
average internal field .B SI1 = (B) — B app at several val- 
ues of -B app are shown in FigJT] for data measured using 
LTF. Consistent data are obtained between runs on LTF 
and MuSR despite the difference in orientation of -B app , 
indicating that the mosaic sample provides a good ax- 
ial average about the a axis. Both -B rms and |2? s h| fall 
continuously with increasing T, reaching zero above a 
characteristic temperature T 1 . The effect of increasing 
B app is clearly seen as a depression of both T" and the 
size of the SC response measured by B rms and |B s h| at 
low T . In order to extract characteristic T' values, fits 
were made for T < T' to a simple two-fluid T depen- 
dence, i.e. 1 — t 4 , where t — T/T'. The obtained values 
are shown in FigOJ; for different values of H. There is 
an offset of around 0.2 K between values of T" obtained 
from B Yms and from -B s h- Previous transport studies^ 
suggested a vortex liquid phase just below T c ; we find 
that T" obtained from B rms is well described by a 3D 
vortex melting curve and we assign this T' to the melt- 
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FIG. 2. (a) Field dependent SC line width at 0.35 K (solid 
points) . Open points are data taken from Refs. [I0I - IT2I The 
global fit is described in the text (solid lines), (b) Comparison 
with B lms (H) for the typical ft-ET organic superconductor. 
Dashed lines show the corresponding Barford-Gunn^ field- 
independent values, (c) Scaling plot of T c against p(0) for 
molecular SCs. The filled point is the present result, open 
points are previous data for molecular SCs^. 

ing temperature T m . The melting curve at low H takes 
the formal T m (H) = T c (0)/(1 + (H/H ) 1/2 ) where H is 
a characteristic field. In highly 2D systems such as k-ET, 
field-induced layer decoupling transitions dominate over 
3D meltin g 18 ! 22 , but strong Josephson coupling stabilises 
the layers for TMC. 

Having established the stability of a VL phase over a 
range of T and -B app , we focus on the accurate determina- 
tion of A. A reliable A value can be obtained by measuring 
-Brms (-Bapp) at low T and comparing with the prediction 
for an ideal vortex lattice calculated using the Ginzburg- 
Landau (GL) model—. The expected behaviour is an 
increase of -B rms with increasing -B app in the region of 
-B c i, followed by a maximum and then a fall towards 
zero at B c2 . The measured B rms (B app ) at T = 0.35 K is 
shown in Fig. [2^, along with the earlier reported results. 
For fitting the data we assume an angular variation for 
£ and A of the form (cos 2 9 + sin 2 6*/7 2 ) -1 / 4 , where 9 is 
the angle to c* and the anisotropy factor 7 is obtained 
from the B c2 slope ratio& as (2.3/0.11) 1 / 2 = 4.5. The 
overall behaviour follows the GL form, consistent with a 
well defined VL at low T. Data above 1 mT are fitted 
by the parameters A = 0.92(2) fim and k = 21(3). The 
corresponding B c2 value of 0.16(4) T is fully consistent 
with direct measurements^. The results of Refs. fTol and 
IrH which were made on aligned crystals with B app || c* 
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FIG. 3. TF data at 1.9 mT taken down to 20 mK (solid 
circles). The square of p(T) normalised to p(0) is plotted 
along with an empirical power law fit (solid line). The linear 
dependence for states with line nodes on the FS (dot-dash 
line) is not seen in the data, which are much closer to fully 
gapped BCS behaviour (dashed line). The inset illustrates 
the anion-ordered FS in the ab plane calculated by Nagai et 
alpl; our data are consistent with a gap function A(k) having 
a nodal plane parallel to the warped FS sheets (dashed line), 
e.g. A(k) = Aosinfc a . The ZF relaxation rate (open circles) 
shows no significant change below T c . 

are also consistent with our fit and our analysis includes 
these data. In contrast, the data points of Ref. [HI have 
an anomalously large -B rms and are excluded from the 
global fit. The small B C 2 for -B app ||c* leads to a remark- 
ably strong field dependence of B Tms in our measurement 
range and the extended plateau region typical of strong 
type-II (large k) materials is not present here. In con- 
trast, k-ET (Fig. 2b) has a more well-developed plateau 
due to its larger values of k and B C 2, along with an en- 
hanced -B rms scale resulting from its smaller A. 

Breaking of time-reversal symmetry (TRS) in some 
types of unconventional SC can lead to spontaneous in- 
ternal fields^ 7 * in zero applied field (ZF). These weak 
fields have been observed by pSK in several examples, 
e.g. S^RuCf^ and more recently LaNiCsr^. The spon- 
taneous fields lead to an increase in the magnitude of the 
ZF depolarization rate of the muon polarization in the SC 
phase. To search for this effect, we took ZF /iSR data as 
a function of T, scanning through T c . Measured spectra 
were fitted to an exponential relaxation function e~ AzF * 
and Azf is shown in Fig. [3] (open circles) . Within exper- 
imental error we observe no evidence for a spontaneous 
local magnetic field, consistent with an earlier report^. 

In order to gain further information on the SC pair- 
ing symmetry present in TMC we have examined the TF 
broadening down to 20 mK, taking relatively high statis- 
tics LTF data in a field of 1.9 mT. This field was chosen 
to be close to the maximum of B lms and a field at which 
the VL is thermally stable over a wide range of T. The 



square of the superfluid density p(T), normalised to its 
low T value p(Q) , is shown in Figj3] (solid points) , each 
point corresponds to ~10 7 analysed muon decay events. 
The observed reduction of p(T) with increasing T reflects 
the excitation of quasiparticles, which is highly sensitive 
to the quasiparticle density of states (QDOS) and thus 
the presence of any gap nodes on the Fermi surface (FS). 
A good description of the data is provided by the em- 
pirical power law p(T)/p(0) = 1 - t n (t = T/T m ). Our 
fitted value is n = 2.5(3), which is not consistent with 
the n = 1 linear behaviour predicted for models with FS 
line nodes and seen previously, for example, in /iSR stud- 
ies of k-ET— The clear saturation of p(T) at low T 
is strongly indicative of a fully gapped SC with no FS 
nodes. This fully gapped state agrees with a previous 
thermal conductivity study 3 - 1 -, but differs from the con- 
clusions of a recent ab plane angle-resolved heat capacity 
study^, in which FS gap nodes were inferred^^ from dip 
structure found in C/T data at 0.14 K and 0.3 T, when 
field was oriented at ±10° to a. The only nodal symme- 
try that would be consistent with our data, however, is 
a p x state having no nodes crossing the FS (Fig. 3 inset). 

From our p(T) dependence we can extrapolate to 
T = 0, giving A(0) = 0.86(2) ^m. Scaling behaviour 
T c oc p m has been explored for many classes of SC, in- 
cluding cuprate a 33 ' 34 and pnictides 3 -^, where m = 1, and 
molecular systems 3 -^ where m = 3/2. The new p value 
places TMC significantly below the molecular SC trend 
line 3 ^ (Fig. Et). Introducing non-magnetic impurities to 
TMC strongly suppresses X ^ 37 ' 38 , but regular samples are 
extremely pure. The low T c therefore suggests a different 
pairing mechanism from the other molecular SCs. 

In order to systematically assign the SC phases, we 
summarize in Table 1 key experimental properties, both 
in the low field L-SC state studied here and also in the 
high field H-SC state established for in-plane fields above 
1.5 T (Table 1, upper section), along with the eight 
specific theoretical possibilities for a P-l triclinic sys- 
tem with inversion symmetry, previously identified in the 
group theory analysis of Powell 3 ^ (Table 1, lower sec- 
tion). For the L-SC phase, it can be seen that the only 
assignment consistent with experiment is the OSO state 
(odd- frequency singlet with odd-parity). 

Although early theoretical studies suggested that such 
an odd-frequency SC would be intrinsically gapless 15 , 
subsequent work has indicated that gapless behaviour 
is not an essential feature. In particular, a model 
based on a condensate of composite bosons made up 
of a Cooper pair and a magnon has recently been ex- 
plored in some detail 4 ^. Depending on model parame- 
ters, the QDOS for this model could range from the stan- 
dard sharply-peaked form of the fully-gapped Bardeen- 
Cooper-Schrieffer (BCS) model through to a form with 
a completely closed gap. In an intermediate regime the 
QDOS peak becomes broadened and weakened compared 
to BCS. As p(T) falls faster than BCS in the region above 
0.2 K (Fig. 3), this difference would be consistent with the 
broadening of the QDOS into the gap region expected for 
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TABLE I. Properties of the superconducting states of TMC, experimental (top two lines) and theoretical (based on Powell). 
Bold text marks our main conclusions and unique assignment for L-SC. (E=even, 0=odd, S=singlet, T=triplet, HS=Hebel- 
Slichter, "experimental result or symmetry requirement from theory, b this study, c direction dependent). 



the intermediate regime in this type of model 4 ^, although 
we note that thermal vortex motion could also lead to 
-B rms falling with T in a region just below T m , giving an 
apparent extra fall in p(T). 

Turning to the H-SC phase we see that there are two 
possibilities, either an ETO state (even-frequency triplet 
with odd-parity) of the polar form, as previously sug- 
gested by Powell 39 , or alternatively an OTE state (odd- 
frequency triplet with even-parity), also of polar form. 
For a quasi-lD system such as TMC a field-induced 
transition from OSO to OTE pairing has been found 
in theoretical studies based on the extended Hubbard 
model, due to enhancement of charge fluctuations over 
spin fluctuations^— . Clear discrimination between the 
two H-SC candidates could be gained by studying the ef- 
fect of disorder on T c within the H-SC phase, however the 
studies to date have only measured the L-SC phas e 37 ' 38 . 

Powell 3 ^ has suggested a scenario in which the field- 
induced transition is due to spin-orbit coupling (SOC). 



In this picture a strong-SOC p-BW stated for the L-SC 
phase crosses over to a weak-SOC p-polar stated in the 
H-SC phase. However, significant SOC is not expected 
here at low T due to the underlying inversion symmetry 
once the TMTSF torsional modes have been frozen out4£. 

Finally, we note that the thermodynamic stability of 
odd-frequency bulk SCs was initially questioned, but re- 
cent theoretical work now concludes that these states 
should actually be stabl o 46 ' 47 . The example identified 
here would suggest that such an odd-frequency state can 
indeed exist in a bulk material. 
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